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ABSTRACT 

We describe a simple algorithm for classifying orbits into orbit families. This algorithm 
works by finding patterns in the sign changes of the principal coordinates. Orbits in 
the logarithmic potential are studied as an application; we classify orbits into boxlet 
families, and examine the influence of the core radius on the set of stable orbit families. 
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1 INTRODUCTION 

The dynamical structure of an elliptical galaxy may be rep- 
resented by a collection of orbits (e.g. Schwarzschild 1979). 
In this representation the time-average of each orbit yields 
an invariant density distribution, and a weighted sum of such 
distributions equals the overall mass density of the galaxy. 
Orbits are categorized into different families depending on 
their shapes; certain orbit families become stable or dom- 
inant in certain potentials. The orbits available in a given 
potential may be characterized in several interrelated ways. 
Classifying orbits into families says something about their 
geometry and thus about the contribution they can make 
• to the mass distribution. This paper describes a simple al- 
gorithm for classifying orbits in elliptical potentials, and il- 
lustrates its application to the well-studied logarithmic po- 
tential (see Miralda-Escude & Schwarzschild 1989, hereafter 
MES89). 

The orbital content of triaxial potentials has long been 
of interest to dynamical astronomers. Integrable potentials 
have fairly simple structures, with one family of 'box' or- 
bits and three of 'tubes' (de Zeeuw 1985). In nonintegrable 
potentials the situation is much more complex; chaotic or- 
bits may appear (Goodman & Schwarzschild 1981; Valluri & 
Merritt 1998), and the monolithic tube and box families are 
partly replaced by a plethora of minor orbit families (Binney 
& Spergel 1982; MES89; Lees & Schwarzschild 1992). This 
rich structure is created by a pattern of resonances between 
fundamental orbital frequences. In d dimensions there are 
d such frequencies uii, and an orbit family corresponds to a 
resonance condition of the form 

n-w = 0, (1) 

where n is a vector of integers, and |n| > 0; thus each orbit 
family may be characterized by the corresponding vector n. 
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Several different approaches to orbit classification ex- 
ist. Spectral methods (Binney & Spergel 1982, 1984; Pa- 
paphilippou & Laskar 1996, 1998; Carpintero & Aguilar 
1998; Valluri & Merritt 1998) classify orbits by measuring 
their fundamental frequencies uii. These methods are pow- 
erful and general but also somewhat complex to implement. 
Alternately, orbits may be classified by heuristic methods 
based on known characteristics of specific orbit families. For 
example, tube orbits maintain a definite sense of circula- 
tion about a principal axis (e.g. Barnes 1992). Likewise, the 
'boxlet' orbits found in nonintegrable potentials may be clas- 
sified by measuring extrema along the principal coordinates 
(e.g. Schwarzschild 1993). 

We have developed a classification algorithm that cap- 
tures the geometry of both loop and boxlet orbits in a simple 
unified scheme. This algorithm is easily implemented and 
generalized to find higher-order resonances. In its present 
form the algorithm is best suited for two-dimensional po- 
tentials; while it provides a partial characterization of orbits 
in three-dimensional potentials, it cannot recognise the full 
range of orbit families present in such systems. Nonetheless, 
the generality and simplicity of our method may make it at- 
tractive when a full-blown spectral treatment is unnecessary 
or prohibitively expensive. 

In this paper we describe the algorithm and present 
some sample results. Section 2 develops the theoretical ba- 
sis for the algorithm and discuss some aspects of its im- 
plementation. Section 3 illustrates the use of this algorithm 
in exploring the orbits in the logarithmic potential. Conclu- 
sions and further applications of the algorithm are discussed 
in Section 4. 



2 ALGORITHM 

The key idea of our algorithm is to determine an orbit's fam- 
ily from the pattern of sign changes of the principal coordi- 
nates (Barnes 1998, p. 352). This approach can be placed on 
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a firm foundation by examining sign-change patterns for the 
closed, stable orbits which parent an orbital family. A given 
family actually contains a one-parameter sequence of stable 
closed orbits, each with a different energy. In what follows 
we will use a single one of these closed orbits as represen- 
tative of an entire family. This simplifies our exposition but 
has no other consequences, since all the closed orbits in a 
given family have the same topology and therefore the same 
sign-change pattern. 

2.1 Closed orbits 

In two dimensions, the resonance condition (1) establishes a 
unique relationship between the two fundamental frequen- 
cies of an orbit. We write (1) in the form 

£lo x + rauy = , (2) 

where £ and m are integers, and w x and uj y are the orbital 
frequencies. We assume with no loss of generality (i) that 
the frequencies obey < cu x < u) y , (ii) that m < < £, (iii) 
that m and £ are relatively prime, and (iv) that the orbit 
has a period of 2tx. The resonance condition then implies 
that w x = — m and uj y = £. 

Since u> x /co y is a rational number, an orbit with two 
frequencies cj x and u) y will close on itself after a finite time. 
Such a closed orbit, if stable, can parent an orbital family. 
Although the motion is generally not harmonic, these closed 
orbits are topologically equivalent to Lissajous figures, which 
are parametric curves generated by 

x(t) — cos(oj x t) , y(t) = cos(uj y t + <f>) . (3) 

where <f> = fc7r/2u> x is a phase-shift. For any given o; x and ui y 
there are two values of the parameter k which yield possible 
parenting orbits: k — generates the parents of boxlets, 
while k — 1 generates the parents of 'antiboxlets' (MES89). 

Using the parametric equations (3), we can derive the 
pattern of sign changes of x(t) and y(t) for any distinct 
choice for (£, m) and k. Since cos(6 l ) has two sign changes for 
8 6 [0, 27r), this pattern is periodic with a period of 7r, or half 
the orbital period. We refer to the pattern generated over the 
interval t € [0, 7r) as the semi-pattern; in each semi-pattern 
x(t) has — m sign changes and y{t) has £ sign changes. Before 
presenting a general treatment, we will discuss the semi- 
patterns of some key orbits, organized by their £ + m values. 
These examples also serve to introduce our notation for sign- 
change patterns. 

2.1.1 The case £ + m = 

If £ = —m then since I and m are relatively prime we have 
(£,m) = (1,-1), so u; x = uj y — 1. In this case, the curve 
defined by (3) with k — is trivial: a diagonal line passing 
through the origin. On the other hand, the curve produced 
with k — 1 is a circle, which is topologically equivalent to a 
loop orbit. For t £ [0, ir) the function y(t) changes sign at 
t = 0, while x(t) changes sign at t = tt/2. Abstracting away 
the actual values of t but preserving the order of events, 
this orbit's semi-pattern is represented by the string 'YX', 
where 'Y' marks a sign-change of y(t) and 'X' a sign-change 
of x(t). An orbit which produces a sign-change string of the 
form 'YX' repeated ad infinitum is a loop orbit. 



2.1.2 The case £ + m = 1 

If I = 1 — m and k = then (3) yields 'centrophobic' orbits 
(MES89), which avoid the centre of the potential. The semi- 
pattern for such orbits is a string of the form YX. . . Y', in 
which — m pairs of the symbols 'YX' are followed by a final 
'Y'. For example, the resonance (£,m) = (2, —1) yields the 
semi-pattern 'YXY', and an orbit which produces a sign- 
change string consisting of infinite repetitions of 'YXY' is 
a centrophobic (2, —1) resonance, also known as a 'banana' 
orbit (Binney 1982). Likewise, the resonance (3,-2) yields 
'YXYXY', which identifies a 'fish' orbit. Unlike the loop or- 
bit described above, all semi-patterns of this kind must begin 
and end with the symbol 'Y'. The sign-change string gen- 
erated by an orbit of this kind will therefore contain places 
where the 'Y' appears twice with no intervening 'X'. These 
doubled symbols have a simple physical interpretation: the 
angular momentum of the orbit reverses between the first 
and second Y'. 

In contrast, k = 1 yields 'centrophilic' orbits, which pass 
exactly through the centre of the potential. The simplest 
such orbit is the 'antibanana', produced by (£, m) = (2, —1). 
For this orbit y(t) has zeros at t = and tt/2, while x(t) 
has a zero at t = tt/2. To represent the simultaneous ze- 
ros of x(t) and y(t) at t = tt/2, we surround the symbols 
with square brackets; thus the semi-pattern for this orbit is 
'Y[XY]\ Next is the 'antifish', produced by (£, m) = (3, -2); 
the semi-pattern for this orbit is 'YXY[XY]\ In all cases 
with £ + m — 1, the semi-patterns for the centrophilic orbits 
resemble those of their centrophobic counterparts, except 
that the final symbols are '[XY]', indicating that the orbit 
passes through the origin. 



2.1.3 The case £ + m = 2 

If I = 2 — m then k — yields centrophilic orbits which pass 
through the origin at times just before t = tt/2; thus the 
semi-patterns for such orbits contain a central '[XY]' pair. 
The simplest orbit of this kind with io x /uj y > 1/2 is (£, m) = 
(5,-3), which produces the semi-pattern 'YXY[XY]YXY'. 
The next resonance in this series, (£,m) = (7,-5), yields 
'YXYXY[XY]YXYXY'. In general, each such pattern has 
three parts: the first is an alternating string of (£ — l)/2 
copies of 'Y' and (— m— l)/2 copies of 'X', the second is the 
central '[XY]' pair, and the third is identical to the first. 

For k — 1 we obtain centrophobic orbits. Each semi- 
pattern contains at total of I copies of 'Y' symbols and I — 
2 — —m copies of 'X'; thus within the semi-pattern two 
'Y' symbols must appear next to each other. The simplest 
example is (£, m) = (5, -3), which yields YXYYXYXY'. As 
already noted, 'YY' pairs indicate a reversal of the angular 
momentum; counting the two 'Y' symbols which begin and 
end the semi-period, we find that a centrophobic orbit with 
£ + m — 2 reverses its sense of angular momentum four times 
before closing on itself. 



2.1.4 General treatment 

The coordinate functions x(t) and y(t) have zeros in the 
interval t £ [0, tt) at times 
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Using these expressions, it's straightforward to generate the 
semi-pattern for any given (£, m) and k. Table 1 gives results 
for all resonances satisfying ui x /ujy > 1/2 which yield semi- 
patterns of 20 symbols or less. The resonances are grouped 
by £ + m value; within each group, they are sorted in order 
of increasing w x /aj y ; finally we list the centrophobic version 
of each resonance before its centrophilic counterpart. The 
column labeled 'Key' gives a one-character symbol for each 
orbit family, which we will use in presenting results of orbit 
classification in Section 3. 

In passing, we note a crucial fact: no semi-pattern con- 
sists of multiple copies of any other semi-pattern. This fol- 
lows because £ and m are relatively prime. 

The motion described by (3) is reversible; that is, the 
transformation t — » — t leaves trajectories unchanged. Thus 
it may seem puzzling that some semi-patterns are not palin- 
dromes. The resolution of this puzzle lies in the fact that 
the transformation t — > t + Ai, where At is a constant, also 
leaves trajectories unchanged. Consequently, all cyclic per- 
mutations of a given semi-pattern should be identified with 
each other since they all represent the same trajectory. The 
reverse of a given semi-pattern can always be cyclicly per- 
muted to obtain the original semi-pattern. We conclude that 
the asymmetric appearance of some semi-patterns is an ar- 
tifact of our arbitrary choice for the origin of the time-line. 



2.2 Family resemblances 

The closed orbits in elliptical potentials are distorted ver- 
sions of Lissajous figures, and each such figure yields a 
unique semi-pattern. We now argue that all the orbits in 
an orbital family yield similar sign-change patterns. Recall 
that in two dimensions, each orbit family is parented by a 
stable closed orbit. Stability is the key which enables an or- 
bit to be a parent; informally, it implies that other orbits 
which start near the parenting orbit remain near it at later 
times. An orbital family is thus the set of all orbits which 
remain near a stable parenting orbit. 

For a centrophobic family the resemblance is simple; all 
the orbits in a given family generate the same semi-pattern 
as does their parenting orbit. This seems intuitively clear for 
orbits which wander only a small distance from the path of 
their parent, but is it true for even the most far-flung family 
members? Yes, because no member of a centrophobic family 
passes exactly through the centre. Consider a continuous 
ensemble of family members which - up until now - have all 
generated the same sequence of x and y sign-changes, and 
are therefore at present all within the same (x, y) quadrant of 
the system. If some of these orbits next have sign-changes in 
x while others have sign-changes in y then the swath defined 
by their orbits must pass across the origin. But continuity 
then implies the ensemble includes orbits which pass exactly 
through the centre, which is impossible since the family is 
centrophobic. 

On the other hand, a centrophilic family includes at 
least one member - the parenting orbit - which passes ex- 
actly through the centre every time. Other family members 
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Table 1. Semi-patterns of 20 symbols or less with LO x /cu y > 1/2. 
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come arbitrarily close to the centre, but generally pass to one 
side or the other; the swath defined by an ensemble of such 
orbits has finite width as it crosses the origin. Thus some 
of members of the ensemble have sign-changes in the order 
'XY' and some in the order YX', while only a infinitesimal 
subset pass exactly through the centre and produce '[XY]'. 
The sign-change sequence generated by a typical centrophilic 
family member is an infinite repetition of its parenting or- 
bit's semi-pattern with one embellishment - where the par- 
ent's semi-pattern contains '[XY]' the member's realization 
may contain either 'XY' or YX' (but not 'XX' or YY'). 
Thus, each '[XY]' indicates a place where the order of the 
symbols is interpreted as ambiguous. 
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This interpretation complicates matters, since the semi- 
pattern of a centrophobic resonance will match the semi- 
pattern of its centrophilic counterpart. Nonetheless, an orbit 
can be assigned to its appropriate family with little ambigu- 
ity - given a sufficiently long sequence of sign-change data 
generated by following the orbit's trajectory. Match the or- 
bital data against each semi-pattern in turn, taking care 
to test centrophobic resonances before testing centrophilic 
ones; the first successful match gives the orbit's family. With 
a finite sequence of sign-change data, it's possible that a 
centrophilic resonance will be miss-classified as centropho- 
bic, but the probability of such an error is small if enough 
trajectory data is available. 

2.3 Implementation 

A fairly straightforward computer program serves to imple- 
ment the algorithm outlined above. This program computes 
the trajectory (x(t),y(t)) of an orbit, recording the sign- 
changes of a: and y; once a total of -/V c h ang e sign-changes have 
been recorded, these are matched against a finite set of semi- 
patterns. A successful match identifies the orbit's family. If 
none of the semi-patterns match the orbit's sign-change data 
then the latter are searched for any strictly periodic pattern 
of length < L max ; this catches all centrophobic families with 
I — m < L max . If no such pattern is found then the orbit 
may be a higher-order resonance, a true (regular) box with 
incommensurate o; x and u> y , or a stochastic orbit; our algo- 
rithm can't distinguish between these possibilities. 

The trajectory-following section of the program is built 
around a 4 th -order RK integrator with a time-step always 
proportional to the local dynamical time. While a higher- 
order integrator and local error control would improve the 
accuracy of the computed trajectories, convergence tests in- 
dicate that most orbits can be robustly classified even when 
integrated with energy errors as large as 0.1 per cent; to 
be on the safe side, we limit the peak-to-peak energy varia- 
tion to 0.02 per cent. After each time-step, the new values 
of x and y are compared with their previous values; if one 
has changed sign then this event is recorded, while if both 
have changed sign then both events are recorded, with linear 
interpolation used to determine which occurred first. The al- 
gorithm does not try to detect passages exactly through the 
origin; such passages are vanishingly rare. 

The semi-patterns used to try and match the sign- 
change data are generated from the (£, m) and k values 
stored in a table like Table 1. When matching the orbital 
data to a semi-pattern we demand symbol-by-symbol equal- 
ity throughout unless the semi-pattern contains '[XY]'; at 
that point the orbital data may contain either 'XY' or YX'. 
One final subtlety about the matching routine is worth men- 
tioning: since the phase of the orbital data is indeterminate, 
it is matched against all cyclic permutations of each semi- 
pattern. Even with this added complication, the amount of 
computer time spent matching orbital data to patterns is 
insignificant; only integer arithmetic is required, and most 
possible matches are quickly rejected. 

Finally, the routine which checks the orbital data for 
arbitrary periodic patterns is quite simple. To test for a 
semi-pattern of length L we take the first L symbols as the 
model and match them against the rest of the data. In this 
matching process it's not necessary to consider cyclic per- 



mutations of the semi-pattern since it automatically has the 
same phase as the orbit. If the match succeeds, the £ and m 
values of the parenting resonance are determined by count- 
ing the number of 'Y' and 'X' symbols, respectively, in the 
semi-pattern. If the match fails then the routine next checks 
for semi-patterns of length L + 1, . . . , L max - 

The availability of this final 'safety-net' for arbitrary 
periodic patterns poses an interesting question: do orbit in- 
tegrations ever yield periodic sign-change sequences which 
can't be generated by the scheme in Section 2.1? No evi- 
dence for such orbits has been found. For example, Table 1 
contains all semi-patterns of length L < 20 which can be 
generated by the scheme. Among the many thousands of 
2-D orbits studied in the next section we have found none 
with periodic semi-patterns of length L < 20 which are not 
in Table 1. This success illustrates the power of MES89's 
ansatz linking periodic orbits to Lissajous figures. 



3 APPLICATION 

To demonstrate the use of our algorithm we made a survey 
of orbits in the logarithmic potential, 

^iln^+^ + lP) , (5) 

with ellipticity b < 1 and core radius R c (Richstone 1980; 
Binney & Tremaine 1987; MES89). The goal of this survey 
was to show how the core radius influences the set of orbit 
families. We focused on orbits of total energy 

£ = $(r) + i|r| 2 =0. (6) 

This restriction to E — is largely a matter of convenience. 
In a singular logarithmic potential (R c = 0) an orbit with 
E 7^ can always be rescaled to one with E = 0. If R c ^ 
then rescaling is still possible if the core radius changes along 
with the energy so that E — In R c stays constant. Thus fixing 
E and changing _R C , as we do in Section 3.2, is equivalent to 
fixing R c and changing E, and no new orbital families would 
come to light if we were to vary both parameters instead of 
just one. 

In designing an orbit survey, perhaps the most critical 
choice is the 'start space' which samples the manifold of 
possible orbits. Below we examine several possible choices. 
First, we present results for a two-dimensional start space 
which is guaranteed to include all orbits but awkward to 
use when surveying a large number of potentials. Second, 
we present results for several one-dimensional start spaces 
which between them allow access to all known orbits. 



3.1 2-D start space 

The phase space for a two-dimensional potential has four 
dimensions: two positions, (x,y), and two velocities, (x,y). 
Since all orbits eventually cross the y axis, we stipulate that 
the initial x coordinate is xo = 0. The equation for E above 
then fixes any one of the remaining coordinates in terms of 
the other two. We constructed a grid in yo and yo, and found 
xo using (6). The resulting two-dimensional start space is pa- 
rameterized by the same coordinates often used when con- 
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Figure 1. Orbit families accessible from the 2-D start space. The potential has core radius R c = and ellipticity b = 0.7. The symbol 
plotted at each grid point indicates the classification according to the 'Key' in Table 1. 



structing a 'surface of section' (Henon & Heiles 1964; Con- 
topoulos 1983; Binney & Tremaine 1987, p. 117). 

Fig. 1 shows the result of classifying orbits in the log- 
arithmic potential with core radius R c — and ellipticity 
b = 0.7 using this two-dimensional start space. We empha- 
sise that this diagram is not a surface of section; rather than 
plotting an orbit's position on the (y, y) plane each time it 
intersects x = 0, we plotted the results of classifying sep- 
arate orbits launched from each point on the grid. We use 
the symbols listed in the 'Key' column of Table 1 to in- 
dicate the classifications, with a small dot marking orbits 
which don't belong to any family with I — m < 20. Our 
result may be directly compared to those of other orbit- 
classification schemes (Carpintero & Aguilar 1998, Fig. 16) 
and to the traditional surface of section for this potential 
(MES89, Fig. 1); the agreement is very good. Basically, the 
phase space for this potential is almost entirely occupied 
by a few centrophobic orbit families. Starting at the lower 
left we find a large number of (£,m) = (2, —1) banana or- 
bits. Moving diagonally across the plot, we next encounter 
a strip, broken in the middle, of (5, —3) resonances. Next 
comes a region of (3, —2) fish orbits, followed by a narrow 



strip of (4, —3) pretzel orbits. Finally, the right-hand side of 
the plot is entirely occupied by (1, —1) loop orbits. 

The classifications shown in Fig. 1 were made using 
Change = 1000 sign-changes per orbit. This gives a high 
level of confidence; for example, every fish orbit in this dia- 
gram has executed 200 semi-patterns without a single sign- 
change out of place. Since the orbit families in this particular 
potential all have fairly short semi-patterns, we could have 
obtained results almost as good with less than 100 sign- 
changes per orbit. A potential which supports a wider range 
of resonances provides a more stringent test of the classi- 
fier. We therefore chose a potential with R c = 0.02 and 
b = 0.8; compared to the potential used for Fig. 1, the finite 
core radius helps stabilize centrophilic orbits, while the more 
modest flattening favors high-order resonances. The same 
orbits were independently classified for iV c hango = 125, 250, 
500, and 1000, and the results compared to determine which 
orbit classifications are sensitive to this parameter. Even 
with Af changc = 125, the (£,m) = (2,-1), (5,-3), (3,-2), 
and (4, —3) resonances were accurately classified, but incor- 
rect classifications sometimes resulted for higher resonances 
such as (9,-5), (7,-4), (8,-5), and (7,-5). We observed 
two kinds of classification errors: false positives from non- 
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Figure 2. Orbit families accessible from the stationary start space. The potential has ellipticity b = 0.8. The vertical axis shows R c , the 
core radius, while the horizontal axis shows /3, the angular coordinate of the initial position. 



resonant orbits which happen to lie near stable resonances, 
and mistaken identities due to confusion between centropho- 
bic and centrophilic versions of the same resonance. The out- 
put of the classifier converges nicely as -/V changc is increased; 
taking the results for -/V ch ango = 1000 as fiducial, we find to- 
tal error rates of 4.3, 1.9, and 0.8 per cent, for iV c hangc = 125, 
250, and 500, respectively. 



3.2 1-D start spaces 

In most previous studies of boxlets the initial conditions 
were generated by choosing initial positions ro on the 
equipotential surface <3>(r ) = E and setting the initial veloc- 
ities ro to zero. This 'stationary' start space yields a wide va- 
riety of boxlets, but it excludes orbits like the antifish which 
don't have a stationary point - that is, a point at which 
r| = 0. Alternatively, if the potential is finite at the origin 
one may consider setting the initial positions ro to zero and 
choosing initial velocities satisfying |ro| = (2E — 2$(0)) 1 ^ 2 
(Papaphilippou & Laskar 1996). This 'central' start space 
includes the antifish, but it excludes centrophobic orbits. 
When using the stationary start space we parameterize the 
initial position in terms of the angle f3 between ro and the x 



axis; likewise when using the central start space we param- 
eterize the initial velocity in terms of the angle f3' between 
ro and the x axis. 

Fig. 2 presents a survey of the stationary start space 
for logarithmic potentials with ellipticity b = 0.8 and core 
radii R c in the range 0.16 to 0.0. As in Fig. 1, the result 
of each orbit classification is shown using the symbols listed 
in Table 1. For the larger values of R c boxlets are few and 
far between, with only the (£, m) = (4, —3) and (6, —5) res- 
onances occupying more than a handful of orbits. As R c 
decreases, low order resonances such as (2,-1) and (3,-2) 
appear and grow in importance, eventually dominating the 
start space as R c — » 0. The general trend with core radius 
seen here is entirely consistent with the results presented by 
MES89. 

Although the overall outline of Fig. 2 is unsurprising, 
some of the higher-order resonances present in this diagram 
may be less familiar. For example, there is a narrow band of 
centrophilic (I, m) — (5, —3) orbits extending from the top 
of the diagram and ending near 7? c — 0.02 and f) — 39°. 
This family is distinct from the centrophobic (5, —3) reso- 
nance seen in Fig. 1 and in previous studies (MES89; Lees 
& Schwarzschild 1992). In addition, this centrophilic family 
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Figure 3. Orbit families accessible from the central start space. The potential has ellipticity b = 0.8. The vertical axis shows R c , the 
core radius, while the horizontal axis shows (i 1 , the angular coordinate of the initial velocity. 



is flanked on either side by two centrophobic higher-order 
resonances: the (7, —4) and (8, —5) families. Only the for- 
mer has been noted in the literature (Lees & Schwarzschild 
1992), but the overall pattern seen here might have been 
anticipated since rational numbers, each corresponding to 
a possible resonance, are dense among the reals. The pat- 
tern of resonant orbits in Fig. 2 may be summarized as fol- 
lows: at some starting angle f3 between any two resonances, 
(<i,mi) and (£2,^2), one may look for a resonance of the 
form (£1 + £2, mi + mi). Of course, not all resonances are 
stable, and unstable resonances don't parent orbit families. 
Nonetheless, we expect that much of the space in Fig. 2 oc- 
cupied by unclassified orbits (dots) is actually threaded by 
a fine array of still higher-order families. 

But Fig. 2 can't capture all the stable resonances; it 
only shows orbits which have a stationary point. In Fig. 3, 
constructed using a central start space, we partly remedy 
this deficiency; this plot includes all centrophilic orbits, both 
those that have stationary points and those that do not. 
We find a number of new orbit families. For larger values 
of R c there is a region of (£, m) = (3, —2) orbits (antifish) 
centred on /3' ~ 45°, and a narrow strip of (5,-4) orbits 
near f3' ~ 78° . At intermediate values of 7? c there is a fairly 



wide region of (4, —3) resonances (antipretzels) near /3' ~ 
75°. Finally, the bottom part of the plot contains regions 
of (7, —4) resonances near ~ 18° and (8, —5) resonances 
near ~ 47°. These families, and several others, are not 
represented in Fig. 2. 

Together, Figs. 2 and 3 exhibit an interesting regular- 
ity; as R c changes, a given boxlet vanishes from one plot at 
almost the same point where the corresponding antiboxlet 
appears in the other, and vice versa. For example, the cen- 
trophilic (£,m) = (3,-2) family is confined to R c J> 0.12, 
while the centrophobic version of the same family is con- 
fined to R c 0.12. A similar statement applies to the pair 
of (4, —3) resonances and to other pairs of resonances as well. 
This behavior can be understood from MES89's discussion 
of bifurcations in the logarithmic potential; they remark that 
as a function of R c , 'a boxlet and its antiboxlet change sta- 
bility in such a manner that generally one is stable and the 
other is unstable'. 

Does this rule extend to all resonances? The centrophilic 
(£, in) = (5, —3) family is seen in both Figs. 2 and 3 for 
Rc iJ 0.02, but neither plot shows the complementary family. 
This is easily explained; the missing family is centrophobic 
and has no stationary point, so it's automatically excluded 
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Figure 4. Orbit families accessible from the perpendicular start space. The potential has ellipticity b = 0.8. The vertical axis shows R, 
the core radius, while the horizontal axis shows yo, the initial position along the minor axis. 



from both plots. We therefore constructed a 'perpendicular' 
start space by setting xo — yo = and using (6) to fix xo in 
terms of yo- This start space allows access to centrophobic 
orbits which don't have a stationary point - provided that 
each orbit at some point crosses the y axis at an angle of 
90°. Fig. 4 shows the results of classifying orbits from this 
start space for yo in the range 0.004 to 0.36. (We excluded 
the axial orbit obtained with yo = 0, and trimmed off the 
loop orbits which dominate this start space from yo — 0.3 
out to the maximum value, y ma x = 6(1 — Rc) 1 ^ 2 — 0.8.) This 
new plot shows the missing centrophobic (5, —3) family ap- 
pearing as expected for R c ^ 0.02; it's centred on yo ~ 0.18. 
Several other families in Fig. 4 also reinforce this pattern; 
for example, the small clutch of centrophobic (11, —7) orbits 
centred on yo — 0.2, R c ~ 0.02 neatly complements the nar- 
row bands of centrophilic (11,-7) orbits which extend down 
to R c ~ 0.03 in Figs. 2 and 3. 

Fig. 4 also shows several families previously undetected. 
The centrophobic (£, m) = (7, —5) resonance occupies a strip 
extending from yo — 0.22 at the bottom of the plot to 
yo ~ 0.26 at the top. This strip has appreciable width for 
smaller core radii but narrows to less than the space between 
adjacent grid points in yo for R c <; 0.08, thereby creating 



apparent gaps where the strip misses the grid. A similar ef- 
fect is seen in the centrophobic (9, —7) family. Both of these 
families seem stable for almost all - if not quite all - values 
of R c tested, so it's probably not worth looking for their 
centrophilic counterparts. 

4 DISCUSSION 

Our algorithm works by reducing orbital trajectories to se- 
quences of symbols, and matching these sequences against 
a set of known patterns. The idea of describing trajectories 
symbolically is not new; it has been applied to geodesies 
on negatively-curved surfaces, to forced van der Pol oscilla- 
tors, and to the restricted three-body problem (Moser 1973 
and references therein). Symbolic representation of orbits 
may have other applications in galactic dynamics. Below we 
briefly discuss a few related issues and some possible future 
developments. 

4.1 Start spaces 

The reader may feel somewhat reassured to learn that the 
three 1-D start spaces considered above are probably both 
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necessary and sufficient to find all boxlet families. We can 
show this by using the properties of the parametric curves 
(3). First, the orbits which have stationary points are those 
for which k = 0. Second, centrophilic orbits are those for 
which the quantity I + m + k is an even number. Third, 
orbits which cross the y axis at an angle of 90° are those for 
which I + 2m + k is even. Now if an orbit has no stationary 
point and is centrophobic then it has k — 1 and I + m + k 
is odd, and because £ and m are relatively prime it follows 
that m is odd and I + 2m + k is even - so the orbit crosses 
the y axis at 90°. Thus any orbit must fall into at least 
one of these three categories. Orbits with k — can be 
accessed from the stationary start space, orbits with even 
I + m + k are reachable using the central start space, and 
orbits with even I + 2m + k result from the perpendicular 
start space. With a few exceptions due to misclassification of 
high-order resonances, these rules accurately predict which 
orbit families appear in each of Figs. 2, 3, and 4. 

The stationary and perpendicular start spaces suffice 
for the singular logarithmic potential (_R C = 0) since all cen- 
trophilic orbits are unstable in this potential (e.g. Kuijken 
1993). 

These results have some bearing on orbit surveys in 
three-dimensional potentials. At a given energy, the mani- 
fold of possible orbits can be covered by a 4-D start space 
constructed by analogy with the 2-D start space described 
in Section 3.1 (Levison & Richstone 1987). Since this start 
space is awkward to sample and to display, many workers fa- 
vor 2-D start spaces analogous to the stationary and perpen- 
dicular start spaces described above (Schwarzschild 1993). 
But these 2-D start spaces don't allow access to centrophilic 
orbits which lack a stationary point, and these orbits may 
be stable in nonsingular potentials. Thus a full survey of 
orbits in a three-dimensional nonsingular potential requires 
at least three different 2-D start spaces: the stationary and 
perpendicular start spaces must be supplemented by a cen- 
tral start space (Merritt 1999). It still remains to be shown 
that these 2-D spaces really are sufficient to access all orbits. 

4.2 Noninteger fe? 

We've assumed, following MES89, that the only allowed val- 
ues of the phase parameter are k = and k = 1. This as- 
sumption is critically important for the discussion in the 
previous subsection. However, D. Merritt (private commu- 
nication) has pointed out that k may have values between 
and 1; for example, in a 2-D harmonic oscillator the x 
and y motions are independent and neutrally stable orbits 
with any real value of k exist. For nonintegrable potentials 
the situation is unclear; k is probably limited to a discrete 
set of values, but it's possible that values between and 1 
are permitted. Our algorithm can't distinguish a resonant 
orbits with noninteger k from a centrophobic orbit with the 
same (£,m), since both yield the same sign-change pattern. 
Such orbits would be excluded from the 1-D start spaces 
used in Section 3.2, though in principle they are accessible 
from a 2-D start space. Visual inspection may be the easiest 
way to discover orbits with noninteger k; such orbits, while 
following the same sign-change patterns as 'normal' centro- 
phobes, would have different symmetry properties. No such 
orbits have been seen in a cursory inspection of candidates, 
but more work is needed. 



4.3 Orbit families 

The results presented in Section 3 emphasize the remarkable 
variety and structure of orbit families in the logarithmic po- 
tential. Of the 46 families listed in Table 1, some 29 can 
be found in Figs. 2, 3, and 4. Almost half of these fami- 
lies are centrophilic. While centrophilic orbits have been de- 
scribed before (MES89; Lees & Schwarzschild 1992; Barnes 
1998, Fig. 42), most earlier studies treated these families as 
isolated curiosities. MES89's bifurcation diagrams showed 
that boxlets and antiboxlets exchange parental roles as R c 
changes; we have made the consequences of this role- reversal 
visible by routinely finding the centrophilic counterparts of 
unstable centrophobic orbits. Centrophilic orbits may be 
destabilized by central mass concentrations in real galaxies, 
but their recognition still helps in understanding the mech- 
anisms which determine the orbital content of a potential. 

We note that some families did not appear in any of 
our classification diagrams. Most of these missing families 
are fairly high-order resonances which may be difficult to 
find because they occupy very small regions of the start 
spaces. Two low-order centrophilic families, however, are 
conspicuously absent. One is the orbit generated by the 
(£, m) — (1,-1) resonance with k = 0, and the other is the 
antibanana, generated by the (2, —1) resonance with k = 1. 
The former resonance is stable in some anharmonic poten- 
tials (de Zeeuw & Merritt 1983), while the latter may oc- 
cur in highly flattened logarithmic potentials (Papaphilip- 
pou & Laskar 1996, Fig. 10). We would not expect these 
centrophilic orbits to appear in our plots since their centro- 
phobic counterparts, the loop and banana orbits, do appear 
and are stable for a wide range of parameters. 



4.4 3-D orbits 

It's somewhat disappointing to admit that the technique de- 
scribed here doesn't trivially generalize to three-dimensional 
potentials. The reason for this failure is fairly straightfor- 
ward. In two dimensions the resonance condition (1) fixes a 
unique relation between the fundamental frequencies, but in 
three dimensions it does not; as a result, the parenting or- 
bits of many three-dimensional orbit families are thin mem- 
branes rather than closed curves (Merritt & Valluri 1999). 
Such parenting orbits don't generate unique periodic sign- 
change sequences, so our present algorithm can't recognise 
them. In contrast, spectral methods work quite well in three 
dimensions (Carpintero & Aguilar 1998; Papaphilippou & 
Laskar 1998; Valluri & Merritt 1998). 

Nonetheless, we can partly classify orbits in three- 
dimensional potentials by projecting their trajectories onto 
the principal planes and using our two-dimensional algo- 
rithm to analyse the resulting sign-change sequences. This 
approach recognises resonances in which one component 
of the integer vector n in (1) is zero; such resonances, 
though not the rule, are fairly common in many potentials 
(Papaphilippou & Laskar 1998; Merritt & Valluri 1999). 
Some orbits obey two independent resonance conditions; 
such 'doubly-degenerate' orbits are parented by closed three- 
dimensional curves. But even doubly-degenerate orbits don't 
always yield strictly periodic sign-change sequences. For ex- 
ample, the centrophobic 'banana-pretzel' orbit, combining a 
(2, —1) resonance in the x-z plane with a (4, —3) resonance 
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in the x-y plane, yields a 'ZYXZ[YZ]X[ZY]ZXYZ' semi- 
pattern containing not one but two ambiguously-ordered 
pairs of symbols. This example suggests that more general 
pattern-matching methods may be able classify the crossing- 
sequences generated by arbitrary three-dimensional reso- 
nances. 

4.5 Future developments 

There are several ways in which this algorithm could be im- 
proved or extended. As just mentioned, one line of attack 
is to generalize the pattern-scanning process. It seems fairly 
easy to extend the routine which searches for arbitrary peri- 
odic patterns (Section 3.2) to recognise centrophilic as well 
as centrophobic semi-patterns; this would provide a better 
safety-net for high-order resonances. It's also straightfor- 
ward, in principle, to couple the trajectory-following and 
pattern-matching routines together so that the latter can 
call the former as needed for more data to resolve ambigu- 
ous cases. This would be more efficient and reliable than 
first integrating the orbit for a fixed number of sign-changes 
and then analysing the results. 

Another line of attack is to adapt the algorithm to anal- 
yse N-body simulations. A simplified version has already 
been used to classify orbits in merger remnants (Barnes 
1998; Bendo & Barnes 2000); however, this version just fol- 
lowed individual trajectories in a frozen approximation of 
the N-body potential. A more interesting project would be 
to classify orbits directly from the computed trajectories in a 
simulation. This could be accomplished by using an N-body 
code to generate separate sign-changes sequences for all N 
bodies; these sequences could be quickly scanned by pattern- 
matching algorithms like the one described here. Such tech- 
niques could attach a concise historical record to each body 
in an N-body simulation, nicely complementing the instan- 
taneous snapshots of positions and velocities which currently 
form the basis for most N-body analysis. 
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